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ABSTRACT 


General differential equations governing the dynamics of the cavitation 
bubbles in a liquid were derived. With the assumption of spherical symmetry 
the governing equations were simplified. Closed form solutions were obtained 
for simple cases, and numerical solutions were calculated for complicated ones. 
The Growth and the collapse of the bubble were analyzed, Oscillations of the 
bubbles were studied, and the stability of the cavitation bubbles were investi- 
gated. The results show that the cavitation bubbles are unstable, and the 
oscillation is not sinusoidal. 
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Section I 
INTRODUCTION 

It is a physical property of any substance that no matter how dense it is, 
there always exist numerous sub-microscopic dimension cavities called nucleii. 

It is from these nucleii that the cavitation bubbles in a liquid begin to grow. 
Inside these nucleii, there may be dissolved gases (such as air) and vapor of 
the liquid. In equilibrium, the total pressure inside the nucleus is equal to 
the local pressure of the liquid at the boundary. When the pressure at the 
boundary falls below the vapor pressure of the liquid at the local temperature, 
the liquid will vaporize, and thus the cavitation bubble will grow. As the 
bubble grows larger and larger, it will become unstable and collapse into small 
bubbles. The growing velocity and limit of the size of the bubble depend on the 
properties of the liquid and the parameters of the flow field. If some of the 
properties of the liquid flow field are fluctuating, the growth of the cavita- 
tion bubble may be oscillating. 

The cavitation bubble will affect the dynamics of the liquid and cause 
damage to the pump and propeller operating in the liquid. In this report, the 
equations governing the dynamics of the cavitation bubble were derived; the 
growth and collapse of the bubbles were analyzed; the oscillations of the bubbles 
were studied, and the stability of the cavitation bubbles investigated. 

It was found that the general equations governing the dynamics of the 
cavitation bubble are difficult to solve without simplification. With the 
assumption of spherical symmetry and other simplifications, the equations were 
reduced to tractable forms. Closed form solutions were obtained for the simple 
cases, and numerical solutions were calculated for the complicated ones. Data 
were presented in both tabulated and graphical forms. The results show that the 
cavitation bubbles are unstable, the oscillation is not sinusoidal, and the 
growing velocity is subsonic in general. 
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Section It 


GENERAL FORMULATION A OF CAVITATION BUBBLES 

2.1 ARBITRARY INTERFACE 

Denote the bubble wall by the closed surface 


S (r , t) = S(x,y,z,t) = 0 


( 1 ) 


The surface S divides the whole space into two regions, R' and R. The 
external region R is occupied by a liquid which may have a substantial amount 
of dissolved gas. The internal region R' is occupied by a mixture of the gas 
and vapor of the liquid. We denote flow parameters in R by unprimed symbols, 
and corresponding parameters in R' by primed symbols. Quantities pertaining to 
the liquid or its vapor will be denoted by the subscript 1, and those pertaining 
to the gas by subscript 2. Symbols without any subscript refer to total 
quantities. Thus for example, the total mole number density n is given by 


n = + n 2 


( 2 ) 


Denote the molecular mass by m, mass density by p , and fluid velocity by 
V ; then we have 


P 1 = m l n l ’ P 2 = m ° n 


p = P 1 + P 2 


2 2 


pv - p^ + p 2 v 2 


(3) 

(4) 


Similar relations hold for primed quantities. 


If the thermal diffusion between the gases is neglected, the governing 
equations in R' may be written as 


3p I ’ 

~ + v * (p i'V> = 0 

9 P2' 

yr + v • <p 2 ,v 2'> - 0 


(5) 


( 6 ) 
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( 7 ) 


1^— + V * (p ’V’) = 0 

p' [f~ + (V' • V)V'J = V • T' + p’b 

P ' [ft + ^ ' ' V)U '] + v • h' = Tr (VV ' • t’) + p ' q ' 

where : 


h' = -A ' VT ' + j kT ' [n 1 ' V x ' + n 2 'V 2 ’ - n’V' ] 


n 1 , n 2 » D 12' [V ~ + 




( 8 ) 

( 9 ) 

( 10 ) 
(ID 



In the foregoing set of equations, V is the gradient operator, s is the entropy, 
T the temperature, p the pressure, U the internal energy per unit mass, h the 
heat flux vector, b the body force per unit mass, and q the heat generated per 
unit time and unit volume. The coefficients of thermal conductivity A, shear 
viscosity y , bulk viscosity 0, and mutual diffusion D ±2 will in general depend 
on the temperature and the relative concentration of the gases. 

In R, one will make use of the fact that p 2 << p i.e., the density of 
the dissolved gas is much smaller than the density of the liquid. 
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Thus one has essentially 

■ h - ? (Sr) ' (r) 


*> -> 
V 2 


-* -> / p 2\ 
pV - pV 1 - DpV ( — ) 


and the governing equations become 


3p l 

sT + v ' <“iV ' 0 


where 


and 


3p 2 -+• f p 2 1 

+ V • (P 2 V 1 ) = DV • I pV — 


at 


+ v • ( P v) = o 


av 


p [_ + (v . v)v] = vx + pb 




P tff + (V * V)U] + V * h - Tr(7V • x) + pq 


-> 

h 


= -AVT 


U = U(T,S) (or some equivalent equation of state) 




k . i 

'k 6 j 


e. . = t-(V. . + V. .) 
ij 2 1,2 j,i' 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


(23) 
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At the bubble wall S, the flow quantities in R and R' may be connected by 
applying the conservation laws across S. Let the unit normal at S(r,t) = 0 be 
N; thus 

N = VS/ | VS j (24) 

and its ith covariant component is 

N i = S ,i / l VS l (25) 


Then, on S(r,t) 

= 0, we 

have 





pi<*i • 

vs + f 

) “ Pi' ( 

V 

•v S + ]f) 


(26) 

p 2 (\ • 

,s + S 

) - DpV 1 


• vs = P 2 ' 

( V 2 ' - vs + |f } 

(27) 

p(v • VS 

+ lt>- 

p* ( V' 

• vs 

+ f) 

at 1 


(28) 


p(V.) ( V • VS + |f) - t. j S 
l 3t ' l 

■ VXj ■ 0S ,i <T ■ 5)> 

p ( u + ) ( V . vs + || )- 

j 2 

= p' ( U’ + ( V' • VS + 

- T ,lj v. 's . + h' • VS + a 

Also from the local thermodynamical equilibrium, we have 

T - T» (31) 


,1 - P’V ± ' ( V • VS+^j 
i = 1,2,3 


(29) 


il -> 

x J V.S . + h 


3S_\ 
9 1 1 

as 

9 1 


i ,3 


(V • N) 


VS 


(30) 
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and 


p 2 = P 2 (P2' ,T) 


(32) 


or 


P 2 = a P 2 ' (33) 

One may remark that a is the surface tension coefficient, and 
^ "N = ^ I" , where and are two principal radii of curvature at the 

point of interest on S = 0. An equation of state is needed to relate the 
partial pressure p' 2 to p' 2 and T' . 

The foregoing then is the general set of governing equations for the 
bubble dynamics. 

2.2 SPHERICAL BUBBLES 

When the bubbles possess spherical symmetry, every physical quantity will 
depend only on the radial distance from the center, and velocity components other 
than the radial one will be zero. Thus the governing differential equations 
will be greatly simplified. The equation for the bubble wall will be 

r - R(t) = 0 (!) 

The unit normal on the bubble wall will be 

7 (2) 

Also 



2 

V • N = - 
r 


( 3 ) 



By use of the general formulations in Section 2.1, it can be shown that the 
governing equations for the vapor and gas inside the bubble are: 


3p l' 9 2 p ’ V 

at + 3r ^l^l^ + 7 = 0 


(4) 



9 2p o ' V o ' 

+ Hr ^2 V 2 ' ^ + r 


= 0 


(5) 


30 ' ■ f- (p'V) + ^L- 0 


3t + 3r 


+ V* 

at 9r 


\ !_ 

av' , 2 v' \ 

/ 3r \ 

3r + r I 


9U1 + v .i!Ll 

at 3r 


P' 

3r 

f 

/ 3V' , 2V'\ 

~ P 

( 3r + r ) 

■ 6’ 

|3T 2V f 

\3r r / 

ah' 

2h’ . , , 

3r 

r + p q 


av' 

3r - 


-f 


- -X' 

3T 1 
3r 

- + f kT 'tVY 

+ w 

- n'V'] 


- V ' 

2 


>2 r 

n n ’ 1 

— I* 1 ') 

n 1 'n 2 



n 1 'n 2 ' 12 |_ 

3r \ n ' / 

+ 

n 

Vp’ 


U' = U' (T T ,S' ) 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


The governing equations for the liquid and dissolved gas outside the bubble 
will be: 


3Pl 

9 

(p x v x ) + 

2pV 




at + 

3r 

= 0 

r 



(12) 

3o 2 

V 1 

3p 2 + 1 

, 2V 

D 3 

r 2 3/ p 2\ 


W + 

3r + P 2 ( 

i77" + r / 

2 3r 
r 

L Pr 3r ( p ) 

(13) 
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( 14 ) 


*£. 

3t 

+ h 

■ (pv) 

+ 2£V 
r 

= 0 

3V 

+ v ■ 

3V _ 

. i iL 

+ ± 

3t 

3r 

p 3r 

P 

/ 

3U t 

3U ^ 


/ 3V 

p ( 

3l + 

v ^) 

|= - p 

( 3r 


9V 

3r 


-?) 


. . /3 2 T ,2 31 \ 

+ M7T + r S7) + P ’ 
3r 


U = U(T,S) 

The boundary conditions at the bubble wall r = R (t) will be: 
P 1 (V 1 - R) = p 1 ’(V 1 ' - R) 

P 2 (Vi - R) - DP f^)- p 2 '(V 2 ’ - R) 


p (V - R) = p ' (V ' - R) 

- p’v<v - a, + p- - ^(fr -f~) 

, / 3V’ 2V' \ 2a 


/3V1 + 2VL\ 2a 

' M 3r + r ) R 


(15) 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


(22) 


+ h ' - r k 


T = T' 


(23) 
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( 24 ) 



Equations (21) and (22) may be re-written as 

p + |2-= p* + p'(V - R) (V- - V) + AlL(U _ x.) 


m + m 

_ 

SV' V' \ 

Ur r / 

3 l 

3r r / 


A ~ + h ' 
3r 


VU _ /9V 

r / p' Ur 


4u 1 / 3V* V* v 0 * / 3V ' 
P r 


2V 


3p ' \ 3r 

1 
2 


: ) 


- \ (R - V) 2 + | (R - V') 2 ] 


(25) 


(26) 


where the latent heat per unit mass for the condensation of vapor is 


l = u'+ e --u- : e - 

P P 


(27) 


2.3 SPHERICAL BUBBLES WITH UNIFORM INTERNAL PRESSURE 


Since the diameter of the spherical bubble is small, the internal pressure 
may be assumed to be uniform. In this case, all physical quantities will be 
essentially functions of time only. The velocities of the vapor and gas inside 
the bubble will be 


V - 


V 


v 

dt 


V ' 
2 


3p 2 ’ 


^2_ 

dt 


V' 


r dp 1 
3p ' dt 


(28) 


The heat flux will be given by 


h' 


( el sQLL _ £_ 1 _ \ 

r l 3 dt ~ 3p' dt ) 


(29) 


In order that the assumption of uniform bubble internal pressure be valid, the 
velocity of the bubble wall must be smaller than the speed of sound and the 
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various speeds of transport in the liquid. The following assumptions can be 
made to further simplify the governing equations. 


(a) Bulk viscosity coefficients are small 


(b) Since the liquid is essentially incompressible and the effect of 
viscosity is small, the interaction between compressibility and 
viscosity may be neglected. 


(c) V = V 1 and p = P 1 

(d) At r = R , V = V x = R 


1 _ ^1 <<c 1 _ ^ 2 . 

P x Dt ^ p 2 Dt 


and 


1 5fi. « 1_ 

p Dt p 2 Dt 


(f ) R 2 « jjp- 


Based on the above assumptions, it can be shown that the governing equations 
for the liquid and dissolved gas outside the bubble are 


n + fe (pV) + * 0 


3P 9 9p 9 

+ V 


St 


Sr 


-»[ 


3 2 P 


Sr 


J + 2^2 

2 r Sr 


— + V — 
St Sr 


1_ JLR 

p Sr 


I SU 
\ St 


+ V^)=-p(f + f) + 


SU 

Sr 




Sr 


4y / SV _ V \ 2 
3 l Sr r / 

+ pq 


(30) 

(31) 

(32) 

(33) 


U = U(T,S) 

At r = R(t) , the boundary conditions will be 
V = R 


(34) 


(35) 
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D 


1 d 


3p, 

Z 

3r 


4ttR 


2 dt 


( 1 " r3 >2) 


(36) 


P 


A 




» + in/ iV V\ 

3 \ 3r ' r / 


( \ TT R 3 p ') 

‘ l+ ^L/ 1V_ Vjl 

\ 3 / 

3p l 3 r r / 


+ £ldu' + 

3 dt p R 


(37) 


(38) 


T = T’ 


(39) 


and 


P 2 ap 2 ' (40) 

The equations governing the uniform vapor and gas inside the bubble will be 
unimportant in the solution of the problem. 
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Section til 

GROWTH AND COLLAPSE OF CAVITATION BUBBLES, INERTIA EFFECT DOMINATING 

3.1 DIFFERENTIAL EQUATIONS GOVERNING GROWTH AND COLLAPSE OF SPHERICAL BUBBLES, 
INERTIAL EFFECT DOMINATING 

When the inertial effect is the dominating factor for the bubble dynamics, 
as in the case of rapid collapse of vapor bubbles in a liquid, the thermal and 
diffusive effects can be neglected in the analysis. In this case, it can be 
shown that the governing differential equations are as follows: 


3p 

3t 


+ 


£ ‘» v > + 


2pV 

r 


= 0 


( 1 ) 


3V , „ 3V 1 3p 

3t + V 3r “ ' p 3r 

At the bubble wall, the boundary conditions are: 


( 2 ) 


V = R at r = R(t) 


and 


P 



= p> + 


4p_ / 3V 
3 \ 3r 



(3) 


(4) 


In general, for most liquids of interest, the pressure will be function 
of the density only: 


P = p(p) 


(5) 


It can be shown that the following solutions satisfy equation (2) 


V = 


2 * 
R R 


( 6 ) 


2 - • 2 

, , , , / R R + 2RR 

p ■ P„V) + P ( — 


rV 

2r 4 


(7) 


3-1 


where p^t) is the pressure field at infinity. The equation governing the 
motion of the bubble wall is determined from the boundary condition, equation 
(4) , to be 


3*2 1 r 

RR + 2 r =p0’ 


(R) 




4pR _ 2q 

R “ R 


( 8 ) 


3.2 SOLUTION OF EQUATION 8, CASE 1 

If the right-hand side of equation (8) is not a function of time, the 
equation can be easily integrated. Let 


P'( R ) “ P OT (t) 


4pR 2 ct 

R " R 


(9) 


then equation (8) becomes 


3*2 

RR + ~ R z 


P 

P 


or 




d_ 
3p dt 


(R 3 ) 


( 10 ) 


Integration from t = t Q to t = t yields 


3*2 

R R 


t 

o 


f R 3 

3p 


t 

t 

o 


Substitution of the limits and simplification gives 



where R = R(t ) 
o o 


( 11 ) 


3-2 



Equation (11) is the simplest solution of equation (8). It shows that the 
velocity of the bubble wall is a function of P/p and the radius R. The graphs 
of equation (11) with P/p as parameter are shown in Figure 3—1. 


3.3 SOLUTION OF EQUATION 8, CASE 2 

If P in equation (9) is assumed to be a function of R by the following 
relation 


P 

P 




( 12 ) 


where R is the radius at t = t , and P some reference pressure. Substitution 
o o o 

of equation (12) into equation (10) gives the differential equation of the 
velocity of the bubble wall as 


RR + 


R 



or 


d 

d t 


3*2 
(R R ) 


2P R ' 
o o 


R 

R 


(13) 


Integration from t 


t to t = t gives 


3*2 
R R 


t 

o 


2P R 
o o 

P 


£n R 


t 

t 

o 


Substitution of the limits and simplification gives 



o 


Equation (14) shows the velocity of the bubble wall as a function of R /R. 

. o 

In equation (14) differentiation of R with respect to R/R q and setting the 

derivative equal to zero gives the value of R/R q for the maximum velocity of 

the bubble wall velocity as: 
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/ 



R_ 

R 


,1 R p 

- exp (7 -2§— ) 


( 15 ) 


F° r R q 10 m/sec, P Q /p = 600 (m/sec) , equation (15) gives 


ir " 1 - 28 

o 

3.4 SOLUTION OF EQUATION 8, CASE 3 

If the coefficient of shear vicosity p is small and neglected, equation 
(8) then becomes 


sS + 2* 2 -p[p’ (R) -r] <17) 

The pressure inside the bubble may be assumed to be uniform and follow Boyle 
Law. 



where p ' = p ' (R ) . 

o o 


(18) 


Also the pressure at infinity p^ is assumed to be constant, 
of equation (18) into (17) gives 



or 


Substitution 


(19) 


, . [2 p ' R 3 ’ „ . 

(R 3 R 2 ) = - ° - ° — - - 2 p R 2 R - 4oRR 

dt p R 


Integration from t = t to t = t gives 


R 3 R 2 


1 

P 


2p ' R £nR - 2p 
00 1 00 3 


- 4a 


2 

Jt 


Substitution of 


the limits and simplification gives 
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( 20 ) 



1 

P 



3.5 SOLUTION OF EQUATION 8, CASE 4 

Here a more general closed form solution will be obtained. In this solu- 
tion, the liquid pressure (p^) and the surface tension (a) are assumed to be 
constant, the shear viscosity is assumed to be small and can be neglected, and 
the gas pressure inside the bubble is assumed to vary as 




n 


or 


P’ 




where p’ = gas pressure at 



(21a) 


(21b) 


R = radius of bubble at t = t . 
o o 

n = constant. 


Substitution of equation (21b) into equation (8) gives 


1 d 

3*2 1 

(RR = - 
P 


R ) 
o' 

n 

1 

2a 1 

2R 2 R dt 

P 0 \ 

R > 


Equation (22) may 

be written 

as 




It cr3r2) 

- ” ° 2R n 
P O 

R 2 ” n 

R 

2p 

1 CO 

p 

2* i 

R R - - 
I 


2p ' 

= ° R n 

P o 

1 

d 

( R 3 " n 

\ 


3-n 

dt 

' 3p 


( 22 ) 
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Integration from t = t to t 

o 


= t gives 


3* 2 
R R 


t 
t 

o 

n 3 


2p ' 0 3-n 

E | 

p o 3-n 


2p “ w 3 2a 2 1Z 

R - — R 

3p p 


( 24 ) 


Substitution of the limits yields 

2p ' 


3*2 3* 2 

R R - R R 


R 3 n - R 3 “ n . 2p 


o o p 


2 R n ( b ( R 3 _ R 3) . 2a 2 . R 2 

o \ 3-n / 3p o p o 


It is convenient to assume that the bubble begins to grow at t = t then 

o 


R = 0 
o 


(26) 


Using equation (26) one obtains 

o f t-» 3-n « 

,3-2 P 2R r W " o ‘ 2p ' 


R R = 


R 

P o 




3-n 


-)- 


- 3 3 

(R J - R ) 
3p o 


f ^ ' R o 2) (27) 


or 


R = 


1 2 p ’o 

r n 

(i) - 

(^) 3 1 

2p 

r oo 

\l -(^) 3 j 

i 2a r 

/ R \ 2 1 

i - — 

) 

| 3-n p 

[} R / 

' R / J 

' 3p 

L ' R / J 

1 " pR L 

' R / J 

1 


1/2 


n 1 3 


For most liquids, the surface tension will be small and may be neglected. Then 
equation (28) reduces to 




(29) 


n + 3 


(25) 


(28) 
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R 


equation (29) 

becomes 


*//P' ^ r R 

, R . o / P 

r R \ 

_ 

/_o\ 1 1Z\ 


▼ ' P ' L R 

' R ' J 3 ' p / 

L 1 ~ l R ) J 


(30) 


When n = 1/2 equation (29) becomes 



(31) 


3.6 THE TIME OF GROWTH 

The time for the spherical bubble to grow from radius (R ) to radius (R) 

o 

can be obtained as follows: Equation (29) can be written as 



(32) 


Integration gives 



(33) 


This integration on the right hand side of equation (33) cannot be carried out 
in closed form and has to be performed by numerical integration. 

3.7 NUMERICAL DATA 

Computer programs were developed to calculate the growing velocity of the 
bubble radius for different values of n and other parameters in equation (29). 
Computer programs were also written to calculate the time of growth as function 
of the ratio of radius to initial radius in equation (33). 
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For 


n = 1 

p /p = 200 (in/sec)^ 

L CO ' 

R =0.02 inch = 0.0508 cm 
o 

the growing velocity R for p'^/p^ = 5, 10, 15 as function of R/R^ are given in 
Table 3-2 and plotted in Figure 3-2; the growing time versus R/R^ ate shown in 
Table 3-3 and plotted in Figure 3-3. 

For 

n = 1/2 

p^/p = 200 (m/sec) ^ 

R =0.02 inch = 0.0508 cm 
o 

the growing velocity R for p' /p = 5, 10, 15 as function of R/R are given in 

o 00 o 

Table 3-4 and plotted in Figure 3-4; the growing time versus R/R are shown in 

o 

Table 3-5 and plotted in Figure 3-5. 
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Growing Velocity R ( m/sec 



Radius Ratio (R/R 0 ) 

FIG. 3.1 BUBBLE WALL VELOCITY VERSUS RADIUS RATIO 
, EQUATION (11) 
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.06 



Radius Ratio (r/R d ) 


FIG. 3.3 GROWING TIME VERSUS RADIUS RATIO 
, EQUATION (33) , n = 1 
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-ABLE 3.1 GROWING VELOCITY OF BUBBLE WALL VERSUS RADIUS RATIO, EQUATION (11) 


Radius n. 
Ratio 

3 0 0 (m/ sec ) 2 

5 00 (m/ sec ) 2 

700 (m/ sec ) 2 

9 0 0 (m/ sec ) 2 

1100 (m/sec) 2 

1300(m/sec) 2 

R/R. 

R(m/ sec ) 

R (m/ sec ) 

R(m/ sec ) 

R(m/ sec ) 

R(m/ sec ) 

R (m/ sec ) 

1 

10 

10 

10 

10 . 

10 

10 

l.ZZ 

12.2 

14. 624 

16 . 701 

18 . 55 

20.225 

21.775 

1.5 

13 .05 

16.254 

18 . 92 

21.26 

23.36 

25.29 

2 

13 . 7 

17 . 44 

20 . 514 

23.184 

25.58 

27 .76 

3 

14.01 

18 . 02 

21.286 

24.114 

26.64 

28.953 

4 

14. 087 

18.16 

21.47 

24.335 

26.9 

29.235 

5 

14.114 

18.21 

21.53 

24 . 413 

26.99 

29.335 

6 

14.126 

18 . 23 

21.56 

24 . 448 

27 . 026 

29.38 

7 

14.132 

18.24 

21.58 

24.465 

27.046 

29.4 

8 

14.135 

18 . 245 

21.586 

24.475 

27.057 

29.414 

9 

14.137 

18 . 248 

21.59 

24.481 

27 . 064 

29.421. 

10 

14 . 139 

18 .251 

21.594 

24.485 

27.068 

29 . 426 

15 

14 . 141 

18.256 

21.6 

24.492 

27 . 077 

29.435 


TABLE 3.2 GROWING VELOCITY VERSUS RADIUS RATIO, EQUATION (29), n = 1 


I 

£ 22 . = 200 (m/s ecf 
P 


> R 0 = 0. 

02 inch = 

D . 0 5 0 8 cm 

Po = c 
Pco 

El = 10 

. Poo ....... 

Bi = 15 

R/R 

R (m/ sec ) 

R/R 

R (m/ sec ) 

R/R 

R (m/sec) 



1.0 

0.00 

1.0 

0.00 

1.0 

0 . 00 

1.1 

11.16 

1.1 

16 . 8 

1.1 

20.98 

1.2 

14.09 

1.3 

23 . 57 

1.2 

26.6 

1 . 4 

16 .28 

1.5 

25.43 

1.4 

31.06 

. 1.6 

16 . 74 

1.7 

25.75 

1.6 

32.28 

1.8 

16.54 

1.9 

25.44 

1.7 

32.37 

2.0 

• 

16.07 

3 . 0 

21.55 

1.9 

3 2.06 

3 . 0 

12.96 

4.0 

18 . 37 

3 . 0 

27.58 

3 . 5 

11.5 

5.0 

15.87 

7 . 0 

16.94 

4 . 0 

10 . 16 

9 . 0 

9 . 29 

11.0 

11.71 

5 . 0 

7 . 73 

11.0 

6 . 86 

15 . 0 

8.11 

6 . 0 

5 . 41 

13 . 0 

4 . 43 

19 . 0 

4 . 91 

7 . 0 

2.65 

14 . 9 

0 .58 

21. 0 

3 . 04 

7.33 

0 . 1089 

14.935 

0.1408 

22.455 

0.1182 

7.38b 

0.00 

14.94 

0.00 

22.46 

0.00 




TABLE 3 . 3 GROWING 


TIME VERSUS RADIUS RATIO, EQUATION (33) 


2 0 0 (m/ sec ) 


0.02 inch = 


t„ ( sec ) 


Po 

Poo 

= 10 

R/R o 

t -t 0 (sec) 


0.0508 cm 


15 


-t 0 ( sec ) 


00000000 

1 

. 00000000 

1 

.00000000 

14610374-02 

2 

. 96031892-03 

2 

.76802543-03 

28431033-02 

3 

.18174659-02 

4 

.22230087-02 

45884556-02 

' 4 

. 28248990-02 

6 

.41223507-02 

68446598-02 

5 

.39986311-02 

8 

. 64853525-02 

99183652-02 

6 

.53505446-02 

10 

.93446914-02 

14965961-01 

7 

. 68965938-02 

12 

.12763646-01 

15356275-01 

8 

. 86595678-02 

14 

. 16850987-01 

15780662-01 

9 

.10672592-01 

16 

.21795669-01 

16246411-01 

10 

. 12985089-01 

17 

. 24690312-01 

16767765-01 

11 

.15674313-01 

18 

. 27958857-01 

17369846-01 

12 

. 18872646-01 

20 

.36198124-01 

18105884-01 

13 

.22841034-01 

21 

. 41866745-01 

19144090-01 

14 i 

.28267025-01 

22 

.50423378-01 

20572564-01 

14 .9351 

.40513069-01 

22.455 

.60710941-01 
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TABLE 3.4 GROWING VELOCITY VERSUS RADIUS RATIO , EQUATION (29), n = 1/ 



O 

CNI 

II 

Jo 

2 

0 (m/ sec ) 

R 0 = 0.02 

inch = 0. 

0508 cm 

Eo = 5 

P 00 

1 1 — ' * — — 

El 

£«_ 

= 10 

Ro 

= 15 

R/Rq 

R (m/ sec ) 

R/R 0 

R (m/sec) 

R/R 0 

R (m/sec) 

| l.o 

! 0.00 

j 

1.0 

0.00 

1.0 

0.00 

1.1 ' 

11.3 4 

1.1 

17.04 

1.1 

21.26 

1.3 

i 16.27 

i 

2.0 

28.54 

1.5 

33 . 98 

1.5 

j 

17 .95 

3 . 0 

27.13 

2.0 

35.78 

r- 

• — i 

18 . 56 

10 . 0 

19.27 

3.0 

34.18 

1.9. 

i 

18.7 

30 . 0 

12 . 6 

10.0 

24.97 

j 3.0 

I 

17.43 

50 . 0 

9.64 

20.0 

20 . 08 

1 4.0,. 

1 

16.01 

70.0 

7.61 

30.0 

17.46 

f— 1 
o 

o 

10 . 91 

90 . 0 

5.94 

100 . 0 

10.33 

! 15.0,.:; 

I 

8 . 55 

110.0 

4.38 

200.0 

6.03 

! 25.0 

5.16 

120 . 0 

j 

3 .57 

300.0 

2.29 

i 3 0.0 

i 

3.56 

140.0 

1.38 

320 . 0 

0. 91 

! 3 5.0 

1.37 

143.0 

0.68 

323 . 0 

0 .45 

3 5.99 

0.06475 

143.995 j 

0 . 04371 

323.995 

0 . 03119 

35.995 

0.00 

144 : : : 

0 .00 

3 24 - 

0. 00 
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TABLE 3.5 GROWING TIME VERSUS RADIUS RATIO , EQUATION (33), n =1/. 


200 (m/sec) 


.02 inch=0. 508 cm 


Po = 5 

Pec 


t - t (sec) 
o 


. 00000000 
. 49730848-02 
. 11158395-01 
. 18883428-01 
.28257247-01 
.33647884-01 
.39583583-01 
.46150076-01 
. 5347o694-01 
. 61730372-01 
. 71226637-01 
.82493776-01 
.96697221-01 
. 11783461+00 
. 14562532+00 


i 16 



143 .9948 


t - t (sec) 
o 


. 00000000 


. 14724970-01 
.37706516-01 
.66972756-01 
.10268568+00 
.14582738+00 
. 17075959+00 


19843756+00 


22945970+00 


26475213+00 

30587830+00 


35588346+00 


.42260022+00 


.52899947+00 


.58399516+00 


= 15 

P oo 


t - t Q (sec) 


00000000 


.33107955-01 
, 84802873-01 
15064630+00 
23099661+00 
32806218+00 
38415800+00 


44643169+00 


. 51622950 + 00 
i .59563499+00 


68816558+00 


80067248+00 


95077552+00 


12341412+01 


13163100+01 


Section IV 


GROWTH AND COLLAPSE OF CAVITATING BUBBLES, THERMAL EFFECT DOMINATING 

4.1 DIFFERENTIAL EQUATION GOVERNING GROWTH OF SPHERICAL BUBBLES, THERMAL 
EFFECT DOMINATING 

The differential equation governing the growth and collapse of spherical 
bubbles is (see equation (8) in Section 3.1): 

rS+ 1“ 2 -?[p'( R ) - - ¥ - ir] <» 

where 

R = radius of spherical bubble 
• dR 

R = = growing velocity of the radius 

p = mass density of the liquid 

p'(R) = pressure inside the bubble, is a function of R, and hence a 
function of time in general 

P^Ct) = Liquid pressure far from the bubble 

y = coefficient of shear viscosity 

a = surface tension coefficient of the bubble surface. 

When the thermal effect is dominating, the effects of viscosity and surface 
tension may be neglected. Thus one obtains the differential equation governing 
the growth of the spherical bubble when the thermal effect is dominating, as 

Rii +|r 2 -i ( p' - pj (2) 

4.2 SOLUTION OF EQUATION (2) FOR CONSTANT VAPOR PRESSURE 

When the temperature of the liquid is kept constant, the vapor pressure 
will be independent of time. If the liquid pressure is also kept constant, 
then equation (2) can be easily solved to yield the growing velocity of the 
bubble wall as 
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where 


R = Radius of bubble at t = t 
o o 

R = Growing velocity of bubble wall at t = t 
o o 


It is convenient to count the time when the bubble begins to grow. Hence 


t = 0 
o 

R q = ° (4) 


Thus, equation (2) reduces to 



The time duration of growth can be obtained by simple integration: 


t = 


r R 


R 


dR 




(5) 


( 6 ) 


4.3 SOLUTION OF EQUATION (2) FOR SLOWLY VARYING PRESSURE 

When the liquid is under heating, the temperature and the vapor pressure 
will be increasing. In this case, p' will be a function of time, and the 
solution of equation (2) will not be so straightforward. However, when the 
heating is slow, equation (2) may be solved by assuming constant p’ without 
introducing much error. Thus all the formulations in Section 4.2 apply to 
this case. 



4.4 NUMERICAL DATA 

Computer programs were developed to perform numerical analysis of equations 
(5) and (6) . For 

p'-p o 

= 300, 500, 700, 900, 1100, 1300 (m/sec) 

P 

the growing velocity (R) versus the radius ratio (R/R ) is plotted in Figure 
4-1 and tabulated in Table 4-1 for equation (5): the growing time (t) versus 
the radius ratio (R/R q ) is plotted in Figure 4-2 and tabulated in Table 4-2 
for equation (6) . 

It is seen that the growing velocity approaches the asymptotic value 
very fast as the bubble grows. 

Although the bubble growing velocity never decreases, the growth will, 
however, eventually stop. Because the cavitation bubble is always unstable, 
the stopping of growth is due to the collapse of the bubble. 
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Velocity R 




3 
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Radius Ratio (R/R 0 ) 


FIG. 4.1 GROWING VELOCITY VERSUS RADIUS RATIO 


, EQUATION (5) 
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Growing Time 



, EQUATION (6) 


4-5 


TABLE 4.1 GROWING VELOCITY 



VERSUS RADIUS RATIO , EQUATION (5) 


700 (m/sec ) 

2 

900 (m/ sec) 

1100 (m/ sec ) 

7 

1300 (m/sec ) 

R(m/sec ) 

R(m/ sec ) 

R (m/ sec ) 

R(m/ sec ) 

0 

0 

0 

0 

15.46 

17 .53 

19.38 

21.07 

CO 

rH 

• 

00 
i — 1 

20.55 

22.72 

24.70 

20.21 

22.91 

25.33 

27 . 54 

21.20 

24.04 

26.57 

28.89 

21.43 

24.30 

26.87 

29 . 21 

21.52 

24 .40 

26.97 

29.32 

21.55 

24.44 

27 . 02 

29.37 

21.57 

24 . 46 

27 . 04 

29.40 

21.581 

24.47 

27.054 

29.41 

21.588 

24.48 

27.062 

29.419 

21.592 

24.483 

27.067 

29.424 

21.599 

24 .49 

27.076 

29.43 










TABLE 4.2 GROWING TIME VERSUS RADIUS 


\ P’-Poo 

RADIUS 

300 (m/sec) 2 

500 (m/sec) 2 

700 (m/sec) 2 

R/R 0 

t (sec) 

t (sec) 

t (sec) 

| .i 

.00000000 

.00000000 

.oooooooo 

2 ' 

.20979336-02 

, 16270816-02 

. 1 3789248-0? 

3 

. 35646598-02 

,27630982-02 

.2339034 1 -02 

4 

.49964936-02 

.38721569-02 

. 32763598-02 

5 

,64187592-02 

.49738390-0? 

.42078511-0? 

6 

.78373242-02 

.60726548-02 

.51361 199-02 

1 

7 

.92541549-02 

. 7 1 701272-0? 

.60636534-02 

1 8 

. 10670065-01 

. 82668862-0? 

.69905840-02 

■ 9 

. 1 20*5840-01 

.93632312-02 

. 79 ! 7 1 647-0? 

10 

. 1 35004*4-01 

. 10859320-0 l 

.88435286-02 

15 

.20572733-01 

. 15937846-01 

. 13473728-01 


, EQUATION (6) 


100 (m/sec) 2 

t (sec) 

• COOOOoOO 

. ! 2! 70827-02 

. 20638894-02 

* 28905 1 9 I -02 
.37! 16645-0? 
.45306735-02 
,534*681 1-02 
. 61 ^ 61570-02 
.69833245-02 

• 7 8003007-02 



. oooooooc 

♦1101 6646-02 
. 18675940-02 
.26153849-02 
.33580988-02 
• 40989203-02 
. 48388360-02 
.557*2708-02 
.63! 78265-02 
.70564093-02 


1300 (m/ sec ) 

t ( sec ) 

.00000000 
.10! 40386-02 
j « 17! *591 1 -02 
.24064219-02 
.30896561 -02 
'.3771 I 128-02 
.445 1 7362-02 
.513 19173-02 
.58 ! ! 841 7-02 
.68916069-0? 


.11883663-01 .10749987-01 .98892031-02 


















Section V 


OSCILLATION OF CAVITATION BUBBLES 


5.1 DIFFERENTIAL EQUATION GOVERNING THE OSCILLATION OF SPHERICAL BUBBLES 

The differential equation governing the motion of the bubble wall of a 
spherical bubble is (see eq. 8 in Section 3.1) 


3-2 

RR + ~ R 



P 


4yR _ 2 g_ 
R R 


( 1 ) 


where 


R = radius of the spherical bubble. 

* dR 

R = = growing velocity of the radius, 

p = mass density of the liquid. 

p'(R) = pressure inside the bubble, is a function of R. 

P^Ct) = liquid pressure far from the bubble, 
y = coefficient of shear viscosity. 

a = surface tension coefficient of the bubble surface. 

Equation (1) is also the nonlinear differential equation governing the 
oscillation of the spherical bubble. 


5.2 FREE OSCILLATION OF SPHERICAL BUBBLES 

Setting the forcing pressure P OT (t) equal to a constant, one obtains the 
differential equation for free oscillation of spherical bubbles as: 


3*2 

RR + j R 


-H’' 


(R) - 


_ _ 

R 


2 a 
R 


( 2 ) 


Since 


3*2 

RR + Tj - R = 


2R 2 R 


d_ 

dt 


3- 2 
(R R ) 


(3) 
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equation (2) may be written as 


^ <*¥> = § (V=) ^(r 2 ,- 


— 4- (R 2 ) - St E R 2 

put p 


( 4 ) 


If the viscosity is small and the last term in eq. (4) may be neglected, one 
obtains 


ft ® 3 * 2 > - ! ( 


p -p 0 


/dt 


< r3 > - r ft < r2 > 


(5) 


Integration from t = t to t = t yields 


3*2 
R R 


JL 

3p 


d 3 2 P 

t P' dT (R )dt “ 3^ 

o 


“ R 3 


— R 2 
P 


3*2 3-2 2 

ft 

, , d , 3 V 

2p 
* 00 

3 3 

R R - R R = — 
° ° 3p 

p'(R) ^ (R )dt 

t* 

" 3p~ 

(R - R o> 

o 

- f CR 2 - R 2 ) 




where 


( 6 ) 


R = radius of bubble at t 
o 

• dR 

R = — at t = t_ 
o dt 0 


t 

o 


Let the time t = t be the time when the rate of increase of the bubble radius 
is equal to zero, namely 


R 

o 


dR 

dt 


t=t 

o 


0 


( 7 ) 


R q = maximum radius of bubble during oscillation. 
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Thus equation (6) reduces to 



where 


Y = ratio of specific heat of the vapor inside the bubble. 

? 

P o = pressure inside the bubble when R = R q 


Substitution of equation (9) into the integral in equation (8) and evaluating 
yields 


_2 

3pR 3 


•t 

P ' (R) (R 3 ) dt = 

' t 

o 


2p : 

(3-3y)p 



( 10 ) 


Substitution of equation (10) into equation (8) gives 

. 1v i ' 

^ - 2P ' 


(3-3y)p 


r/ R v ’ 

/ R V 1 

2P r 

/ R \ T 

— _ 


— Il - 

-2. 

. ' R ' 

' R / J 

3p L 

' R ' J 



(ID 
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5.3 VELOCITY VARIATION DURING OSCILLATION 

Equation (11) gives the square of the growing velocity as: 



where 


( 12 ) 


Rq = Maximum radius of the bubble during oscillation. 

R = Radius of the bubble at any time. 

p = Pressure inside the bubble when R = R 
u o 

p = Mass density of the bubble. 

Pa. = Pressure in the liquid. 

a - Surface tension coefficient of the bubble surface. 

Y = Ratio of specific heat of vapor and/or gases inside the bubble. 

P P 

XT U O 00 O 

Note that p , — and — are non-dimensional and may be designated the 
minimum vapor pressure, the liquid pressure and the surface tension parameter, 
respectively. The surface tension parameter changes as the radius grows. 

For 

P 0 2 

— = 1 (m/sec) 

— = 0.0254 (m^/sec^) 

P 


R Q = 0.02 in = 0.0508 cm 

and 

p 

~ = 50, 100, 150 (m/sec)^ 


(13) 
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the growing velocity as a function of the ratio of radius to maximum radius 
was analyzed. The results are tabulated in Table 5.1 and plotted in Figure 5-1, 
with the growing velocity versus the ratio of radius to maximum radius. It is 
seen that the growing velocity increases very fast and reaches the maximum very 
quickly as the bubble grows from the minimum radius and then decreases very 
fast after the maximum velocity. For the majority of the time, the growing 
velocity is low, as compared to the maximum. At the maximum and the minimum 
radius, the growing velocity is zero. 

As can be seen in Figure 5—1, (or Table 5-1), the maximum growing velocities 
for the three cases are: 

R = 990 m/ sec , 

R = 1456 m/sec , 

R = 1998 m/sec , 


for — = 50 (m/sec) 2 , at = 0.034 
P R o 

Poo 2 R 

for — = 100 (m/sec) , at = 0.028 (14) 

P R 

o 

Poo 2 R 

for — = 150 (m/sec) , at ~~ - 0.024 
p R 


5.4 DETERMINATION OF MINIMUM AND MAXIMUM RADII 

At the minimum and the maximum radius during oscillation, the growing 
velocity vanishes. From Table 5.1, when R = 0 it is seen that minimum radius 
ratios are 


~ = 0.0256 
R 

o 

— = 0.0213 
R 

o 

— = 0.0183 
R 

o 


Poo 2 

for - — = 50 (m/sec) 


for — — = 100 (m/sec)‘ 


for — = 150 (m/sec)‘ 


(15) 


Hence, when the maximum radius is specified, the minimum radius can be 
determined from the minimum radius ratio. Here, the maximum radius is assumed 
to be 


R = 0.02 in = 0.0508 cm 
o 


( 16 ) 
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Hence the minimum radius will be 


= 0.0256 x 0.0508 = 0.00130048 cm , f or = 50 ( m / se c) 2 

R m = 0-0213 x 0.0508 = 0.00108204 cm , for — = 100 (m/sec) (17) 

P 

R m “ 0.0183 x 0.0508 = 0.00092964 cm , for ^ = 150 (m/sec) 2 

5.5 DETERMINATION OF GROWING TIME FROM MAXIMUM TO MINIMUM RADIUS 
To determine the growing time, one writes equation (12) as 

R = F 


or 



(18) 



Integration of equation (19) from the minimum radius (R ) to R = R vield 

m J 

the growing time 



( 20 ) 


fhe integration in equation (20) cannot be carried out in closed form and has 
to be integrated by numerical method. Results of the numerical integration 
of equation (20) for the three cases specified in equation (13) are tabulated 
in Table 5-2 and plotted in Figure 5-2 with growing time versus the ratio of 
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radius to maximum radius. It is seen that the majority of the time is spent 

for — = 0.5 through 1, since the growing velocity for this later part of 
K 

O 

growing is relatively small. 


5.6 DETERMINATION OF PERIOD AND FREQUENCY OF OSCILLATION 

The period of oscillation is equal to twice the growing time from the 
minimum to the maximum radius. It is seen in Table 5-2 that the growing times 
from the minimum to the maximum radius are 


t = 0.001491 sec , for — = 50 (m/sec)' 

P oo 

t = 0.001299 sec , for — — = 100 (m/sec) ‘ 

p» ; 

t = 0.001166 sec , for — = 150 (m/sec) 


( 21 ) 


Hence the periods of oscillation are 


T = 2 x 0.001491 = 0.002982 sec 
T = 2 x 0.001299 = 0.002598 sec 
T - 2 x 0.001166 = 0.002332 sec 


for - — 50 (m/sec) 

P ~ , , 

for — 100 (m/sec) 

P “ , , 

for — 150 (m/sec) 

P 


( 22 ) 


The frequency of oscillation is equal to the reciprocal of the period: 


f -T 


(23) 


Hence the frequencies of oscillation are 


f - 1 

= 335.35 

cps 

•Fr»r* 

t'oo 

— c: n 

(m/sec) 

0.002982 

, ior 

JU 

P 

f - 1 

= 384.91 

cps 

■Fnr- 

p=* 

— i nn 

(m/sec) 

0.002598 

9 ior 

P 

f - 1 

= 428.82 

cps 


Poo 

— I tzri 

(m/sec) 

0.002332 

y ror 

— 1 JU 

p 


( 24 ) 
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Growing Velocity ( m/sec 



FIG. 5.1 GROWING VELOCITY VERSUS RADIUS RATIO, 
EQUATION (12) 
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TABLE 5.1 GROWING VELOCITY VERSUS RADIUS RATIO , EQUATION (12) 



Po 

p 

2 

- 1 (m/ sec ) 

, - = 0.0254 (m 3 / sec 2 ) 

P 

, R q = 0.0508 cm 

2» 

p 

2 

= 50 (m/ sec ) 

P CO 

- o 

= 100(m/sec) 2 

P°o - 

150(m/sec) 2 

R/R 

0 

R (m/sec) 

R/R 

0 

R(m/ sec ) 

R/R 

0 

R(m/ sec ) 

.0256 

. 00000000 

. 0213 

.00000000 

. 0183 

. 00000000 

.0257 

.59016947+02 

.0214 

.16607001+03 

. 0184 

. .19153883+03 

. 03 

.92069067+03 

. 027 

.14457208+04 

.023 

.19763177+04 

.033 

. 98651910+03 

.028 

.14563983+04 

.0 24 

.19980717+04 

. 034 

. 98951041+03 

. 029 

.14541866+04 

.025 

.19949377+04 

.035 

. 98705723+03 

. 04 

.11782545+04 

. 04 

.13816960+04 

.06 

.63052161+03 

.06 

.74286283+03 

.06 

. 84031713 + 03 

.08 

.441538 0'9 + 03 

.08 

.50992818+03 

.08 

.57017332+03 

. 1 

.32795098+03 

' .1 

.37530237+03 

. 1 

. 41731510 + 03 

. 201 

. 12118630+03 

.201 

.13695829+03 

.201 

.15109278+03 

. 302 

.65439157+02 

. 3 02 

.73886077+02 

. 3 02 

. 81461780+02 

. 40 3 

. 41296041+02 

.403 

.46704581+02 

.403 

. 51548742+02 

. 504 

. 28162379 + 02 

.504 

. 31939850 + 02 

.504 

. 35315555 + 02 

. 605 

. 19931228+02 

.605 

.22681422+02 

.605 

. 25132449+02 

• 70 6 

. 14193929 + 02 

. 7 06 

.16212936+02 

.706 

.18006971+02 

• 80 7 

. 97 73 3811 + 01 

. 8 07 

.11207601+02 

.807 

.12478047+02 

CO 

o 

cn 

.58624441+01 

.908 

.67499404+01 

.908 

.75336008+01 

1 

.00000000 

1 

.00000000 

1 

.00000000 1 


CO 
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TABLE 5.2 GROWING TIME VERSUS RADIUS RATIO, EQUATION (20) 




PA- = Km/sec) 2 , £ =0> 02 54 (m 

3 . 2 

/sec ) , R 

= 0.0508 cm 



p 

P 

c 

P oo 

p 

2 

= 50(m/sec) 

p 

100(m/sec) 2 

2a = 
p 

15 0 (m/ sec ) 

R/R 0 

t (sec) 

R/R 0 

t (sec) 

R/R 0 

r~ * 

t (sec) 

.025 6 

. 00000000 

.0213 

. 00000000 

. 0183 

. 00000000 

. 0257 

. 33888829-08 

.0214 

.44848466-08 

. 0184 

. 10441826-08 

. 03 

.16035638-06 

. 027 

.11608848-06 

. 023 

. 70865008-07 

. 033 

.22270341-06 

■ .028 

.12986035-06 

.024 

. 80917655-07 

. 034 

1 .24293612-06 

.029 

.14359555-06 

.025 

. 90927471-07 

. 035 

. 26316590-06 

. 04 

. 30949377-06 

.04 

.26928377-06 

.06 

■ .89061124-06 

.06 

. 74348303-06 

.06 

. 64818597-06 

..08 

.16567892-05 

.08 

.14016946-05 

.08 

.12341341-05 

.1 

. 27164283-05 

.1 

. 23240837-05 

.1 

. 20617139-05 

.201 

.13832176-04 

.2 01 

.12128094-04 

. 2 01 

.10930404-04 

.3 02 

. 37 323278-04 

. 3 02 

. 32931417-04 

.302 

.29797800-04 

.4 03 

.76898401-04 

. 4 03 

. 67954348-04 

.403 

. 61547466-04 

. 5 04 

.13680060-03 

. 5 04 

.12084343-03 

.504 

.10942231-03 

. 6 05 

. 22270192-03 

. 6 05 

.19645095-03 

. 6 05 

.17772590-03 

. 7 06 

. 34332052-03 

.7 06 

. 30229993-03 

. 7 06 

.27313415-03 

.807 

. 51501863-03 

. 8 07 

.45230314-03 

. 8 07 

. 40801848-03 

.90S 

.77959753-03 

.908 

. 68251800-03 

. 9 08 

. 61451709-03 

1 

. 14910834-02 

1 

.12988718-02 

1 

.11659031-02 











Section VI 


STABILITY OF CAVITATION BUBBLES 


6.1 DIFFERENTIAL EQUATION GOVERNING BUBBLE STABILITY 

If some infinitesimally small deviation from the spherical symmetry is 
given to the spherical bubble, one would like to know whether this perturbation 
will grow or diminish. The stability equations are thus the linearized forms 
of the equations that govern the dynamics of the spherical bubble. To analyze 
the stability problem, the perturbation is expressed in terms of spherical 
harnomics. Thus the equation of the bubble wall is: 


r = R + Y a Y 
Ll n n 
n 

where 

R = R(t) = the unperturbed radius of the bubble. 

Y n = spherical harnomic of degree n 

a n = a n (t) = the amplitude of perturbation 


and 



<< R 


( 1 ) 


( 2 ) 


It can be shown that the stability equation for a^ is 


d 2 a 


n 


dt 


R 

R dt 


A a = 0 
n n 


where 


A 

n 


[n(n-l)p - (n+1) (n+2) p * ]R - (n-l)n(n+l) (n+2)q/R 2 
R[np + (n+l)p ' ] 


(3) 


(4) 


It is noted that the mode with n = 1 corresponds to a translation of the 

bubble as a whole, and only the modes with n >_ 2 are relevant to the distortion 

of the spherical bubble. The density of the vapor (p ’ ) , being much smaller than 

the density of the liquid (p) , may be neglected, and the expression for A reduces 

n 
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(n-1) (n+1) (n+2) 


a 


(5) 


A = (n-l)R 



By the substitution of 



t 

o 


where R = the value of R at t = t . 
o o 


Equation (3) can be transformed to 



n 


+ 


B b = 0 
n n 


where 


( 6 ) 


(7) 


B = (n-1) (n+1) (n+2) — — r- - (n - -^ /2)R - |(|f 

n 3 R 4 \ R / 

pR 

6.2 CRITERIA FOR STABILITY OR INSTABILITY 

It can easily be shown that the criteria of stability of differential 
equations in the type of equation (7) are 


E 

j> 0 

stable 


n 

B 

n 

< 0 

unstable 

(9) 


Thus these are the criteria of the stability and unstability of the spherical 
cavitation bubble. Here, instability means that the amplitude of disturbance 
is increasing expotentially , while stability means that the amplitude of dis- 
turbance is decreasing expotentially. 
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6.3 DISCUSSION OF THE CRITERIA 

(a) For n >_ 2, the first term in the expression of B^ in equation (8): 
(n-1) (n+1) (n+2) 

PR 

is always positive. Hence the effect of the surface tension (a) is always 
stabilizing, and this affirms common physical reasoning that the surface 
tension is to prevent the bubble from collapsing. 

(b) When the bubble is growing: 

R > 0 


and when the growing velocity increases with time 

(R) - R > 0 

( 10 ) 
( 11 ) 

Since the surface tension of the vapor-liquid interface is small for many 
liquids; when the bubble is growing, and the growing velocity increases rapidly 
with time, as in the case of boiling, the situation will be unstable. 

6.4 REGIONS OF STABILITY AND INSTABILITY 

For the convenience of parameter study, equations (10) and (11) are 
reduced to non-dimensional forms: 

Stable: 


2 

OT 



the criteria of stability and unstability will be 

. 2 

stable: (n-1) (n+1) (n+2) > (n+|) | + | (|) 

pR 

unstable: (n-1) (n+1) (n+2) < (n+|) ^ ( j|) 

pR 4 ' ' 
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Unstable: 



< F 


where 


F ■ **** ' (— 1) (n+l) (n+2) ~ H) T + I (f 1 ) 2 ]^) 3 

L J o 

R o ~ ra di- us °f bubble at some specific time (for example initial time) 

x - some characteristic time of the bubble (for example the period of 
oscillation) 



with acceleration parameter — - — = 2 

R 

Rt 

velocity parameter — = 10 

R 


degree of harmonic n = 2, 3, 4 


value of function F versus — are given in Table 6.1; the curve of function F 

versus ~ and the region of stability and instability are shown in Figure 6.1. 
o 


6.5 


BUBBLE SIZE AT COLLAPSE 

Equation (28) in Section 3 gives the square of the growing velocity as 


R 2 - 


2 P'o 

- r y 

( _o\ 

( e o) 3 i 

2 P r 

00 -i 

R 

(_o\ I 2o 



(3-y)p 

\ R ' 

’ R ' J 

3p L 1 

\R ) J pR L 1 

\R ) J 

(12) 


where 

p' = pressure inside the bubble 

p^ = pressure in the liquid 

R q = radius of bubble at initial time 

y = ratio of specific heat of vapor gas inside the bubble 
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It can be shown that when the right hand side of equation (12) becomes 
negative, the bubble will collapse. Equating the right hand side of equation 
(12) to zero, one obtains the expression relating the bubble size at the time 
of collapse to other parameters as 


2p '° U- 

T R J 

°\ (_o\ 1 

2 P 

oo 

r, . (h>i 1 + a [ 


(3-y)p L' R 

/ 'R ' J 

3p 

L 1 u ) j + pR l 

1 'r / J 


Equation (13) can be written in non-dimensional form as 




o o 


(14) 


The surface tension term in equation (14) is small and may be neglected. By 
doing so, one obtains the pressure ratio versus the bubble size at collapse as 



(3-y) 

3 



(15) 


For specific heat ratio 
Y = 1.4, 1 

P' Q R 

the curve of the pressure ratio versus the bubble radius ratio — at 

o 

collapse is shown in Table 6.2 and Figure 6.2. 
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Pressure Ratio (p'/p 


it 



AT COLLAPSE 
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TABLE 6.1 



EGIONS OF STABILITY AND INSTABILITY 


« 2 • 


n = 2 


F(R,R,R) 


00000000 

.66666667+01 


.18000000+03 
.42666667+03 
.83333333+03 
. 14400000+04 
.22866667+04 
. 34133333+04 
.48600000+04 
.66666667+04 


o 
1 — 1 

II 

n = 3 

n = 4 

F(R,R,R) 

F(R,R,R) 

00000000 

00000000 

. 20500000 + 01 

. 93333333+00 

.16400000+02 

.74666667+01 

.55350000+02 

.25200000+02 

.13120000+03 

. 59733333+02 

.25625000+03 

. 11666667+03 

.44280000+03 

. 20160000+03 

.70315000+03 

.32013333+03 

.10496000+04 

.47786667+03 

.14944500+04 

.68040000+03 

. 20500000 + 04 

.93333333+03 










Radius 

Ratio 

Y = 1 

Y = 1.4 

Radius 

Ratio 

R/R 

0 

El 

Poo 

El 

P~ 

R/R 

0 

1 

.00000000 

.00000000 

ii 

2 

.15555556+01 

.18377832+01 

12 

3 

.21666667+01 

. 28891200 + 01 

13 

4 

. 28000000+01 

.41027200+01 

14 

5 

. 34444444+01 

. 54508613 + 01 

15 

6 

.40952381+01 

.69152000+01 

16 

7 

.47500000+01 

.84837333+01 

17 

8 

. 54074074+01 

.10147200+02 

18 

9 

.60666667+01 

.11893333+02 

19 

; 


10 


. 67272727 + 01 


.13728165+02 


20 


RATIO (R/R 0 ) AT COLLAPSE 


Y = 1 

Y = 1.4 

El 

El 

Pec 

. Pee 

.73888889+01 

. 15632000+02 

.80512821+01 

. 17610667 + 02 

.87142857+01 

.19653333+02 

. 93777778 + 01 

. 21760000 + 02 

.10041667+02 

. 23940736 + 02 

.10705882+02 

.26133333+02 

.11370370+02 

.28426667+02 

. 12035088+02 

.30773333+02 

. 12700000+02 

.33173333+02 


. 13365079 + 02 


.35645067+02 












Section VII 


CONCLUSIONS AND RECOMMENDATIONS 


7.1 CONCLUSIONS 

The general differential equations governing the dynamics of the cavitation 
bubbles are complicated and difficult to solve. However, with the assumption 
of spherical symmetry and some further simplification, the equations were reduced 
to tractable forms. In this report, the growth and collapse of cavitation 
bubbles were studied, the oscillations of the bubbles were analyzed, and the 
stability criteria investigated. It was found that the cavitation bubbles are 
unstable in general; that the frequency of oscillation is high, and the oscilla- 
tion not sinusoidal; that there is much variation in the growing velocity, and 
that the growing velocity is subsonic in general. 

7 . 2 RECOMMENDATIONS 

It is considered that the theory and results on cavitation bubbles developed 
in this study can be advantageously applied to investigate the effect of 
cavitation bubbles on the dynamics of propellant feed lines in the space rocket 
engines. Since the cavitation compliance of the turbopump is important in the 
analyses of the structural and propulsion systems coupling longitudinal stability 
(POGO), the following are recommended for future study: 

(a) Establish a mathematical model and derive the dynamic equations of 
the propellant feedline system of aerospace rocket engines with the 
effect of cavitation bubbles included. 

(b) Solve the dynamic equations and determine the dynamic characteristics 
of the propellant feedline system. 

(c) Investigate the effect of cavitation bubbles on the dynamics of the 
feedline system and on the POGO problem. 
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Section VIII 


APPENDICES 


8.1 APPENDIX A. NOMENCLATURE 


Notation 

a = a (t) 
n n 


b 


b 

n 



h 

k 

in 


N 


n 


P 

o 

P 

p'(R) 


q 

r 

R 


R = R(t ) 
o o 



R 

o 


Description 

Amplitude of perturbation of bubble radius in nth 
harmonic 

Body force per unit mass 



coefficient of mass diffusion 
Frequencies of Oscillation 
Heat flux Vector 
Boltzmann constant 
Molecular mass 

Unit normal at the bubble wall 
Number density 
Some reference pressure 
Pressure 

Pressure inside the bubble 
Pressure inside the bubble when R = R 

o 

Pressure field at infinity 

Heat generated per unit time and unit volume 

Eulerian spatial radius 

Radius of bubble 

Radius of bubble at t = t 

o 

Growing velocity 
Growing velocity at t = t 
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APPENDICES (Concluded) 


Notation 



s 

T 

t 

U 

V 

Y 

n 

6 

Y 
A 

y 

p 

a 

T ' 

V 


Description 

Growing acceleration 

Minimum radius of bubble 
Entropy per unit mass 
Temperature 
Growing time 

Internal energy per unit mass 
Growing velocity 
Spherical harmonic of degree 
Bulk viscosity 

Ratio of specific heat of the vapor inside the bubble 
Coefficients of thermal conductivity 
Shear viscosity 
Mass density 

Surface tension coefficient 
Component of the tensor 
Gradient operator 
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